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Maximizing Entropy by Minimizing Area: 
Towards a New Principle of Self-Organization 
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We propose a heuristic explanation for the numerous non-close-packed crystal structures observed 
in various colloidal systems. By developing an analogy between soap froths and the soft coronas of 
fuzzy colloids, we provide a geometrical interpretation of the free energy of soft spheres. Within this 
picture, we show that the close-packing rule associated with hard-core interaction and positional 
entropy of particles is frustrated by a minimum-area principle associated with the soft tail and 
internal entropy of the soft coronas. We also discuss these ideas in terms of crystal architecture and 
pair distribution functions and analyze the phase diagram of a model hard-sphere-square-shoulder 
system within the cellular theory. We find that the A15 lattice, known to be area minimizing, is 
favored for a reasonable range of model parameters and so it is among the possible equilibrium states 
for a variety of colloidal systems. We also show that in the case of short-range convex potentials the 
A15 and other non-close-packed lattices coexist over a broad ranges of densities, which could make 
their identification difficult. 



PACS numbers: 61.50.Ah, 82.70.-j 



I. INTRODUCTION 

Colloids are all around us. From milk to microreac- 
tors, from pie filling to paint, and from suspensions to 
sieves, colloidal materials are important in a variety of 
products and technologies and are the basis for a new 
class of functional materials. The advent of applications 
based on engineered crystalline materials, such as pho- 
tonic bandgap materials,E10 has put even greater focus on 
colloids. Both the size of colloidal particles and the inter- 
particle interaction are tunable, which provides the basis 
for the manufacture of ordered structures with desired 
lattice spacings and space groups as well as mechanical, 
thermal, and electrical properties. 

Over the years, considerable amounts of experimental 
crystallography data have been collected for various col- 
loidal systems. LI While the tabulation of the relationship 
between the interparticle interaction and the symmetry 
of the crystal lattice is undoubtedly valuable, it is de- 
sirable, if not necessary, that the synthesis of colloidal 
crystals with specific space group be based on analytic 
rather than empirical insight into the mechanisms of self- 
assembly. In principle this should be possible because the 
interaction between the particles is much simpler than 
in atomic or molecular crystals: colloidal particles are 
made of hundreds or thousands of atoms, and the effec- 
tive interaction between the particles is less specific and 
thus much simpler than interatomic interactions. Thus 
the interaction between colloidal particles is chiefly deter- 
mined by their mechanical as opposed to chemical struc- 
ture. Here we will address those colloidal particles that 
are characterized by a relatively dense core and a fluffy 
corona. In this case, the interparticle potential may be 
approximated fairly well by the simple hard-core repul- 
sion dressed with a repulsive short-range interaction of 
finite strength. 



Present theoretical understanding of the phase be- 
havior and stability of the various colloidal systems is 
quite impressive, especially in the case of charged col- 
loidal suspensions interacting via screened Coulomb po- 
tentials.QLj It js, now well established both by analyti- 
cal approachesu'D and by simulationsQ that the solid part 
of the phase diagram of charged colloids includes face- 
centered cubic (FCC) and body-centered cubic (BCC) 
lattices. In the case of a convex interparticle potential 
- its simplest variant being the square shoulder inter- 
actions - the phase diagram is even more complex and 
includes the dense and the loose FCC and BCC lattices.EI 
More elaborate soft potentials such as the interaction be- 
tween star polymers lead to body-centered orthorhom- 
bic (BCO) and-diamond lattices in addition to FCC and 
BCC lattices.EI 

While extremely important, these theories do not pro- 
vide a robust explanation of the stability of colloidal crys- 
tals, and the aim of this study is to look at the problem 
from a more geometrical point of view and to capture the 
statistical mechanics of colloids in a new self-organization 
principle. One such principle is the maximum packing 
fraction rule, which states that pure excluded-volumc 
interactions favor an expanded close-packed structure, 
thereby maximizing the configurational entropy. In the 
case of monodisperse hard sphexes—such an arrangement 
corresponds to the FCC lattice.LxE3 However, many soft- 
sphere systems form non-close-packed lattices, including 
the BCC, BCO, and diamond lattices as well as the A15 
lattice&J observed in crystals of self-assembled micelles of 
some dendritic polymers O Can we understand the ex- 
istence of this rather loosely packed structure? Is there 
another simple geometrical principle that describes the 
self-organization of soft spheres and is analogous but op- 
posing to the maximum packing fraction rule? We pursue 
this question using an idealized model and find an anal- 
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ogy between the soft-sphere crystals and dry soap froths, 
the latter being described by Kelvin's problem of find- 
ing the minimal-area regular partition of space into cells 
of equal volume. Within this framework, we propose a 
new principle of area-minimization that can favor these 
loosely packed lattices of soft spheres.li3 

Having proposed our semiquantitative and intuitive 
explanation for the behavior of soft-sphere colloids, we 
have checked our predictions within a more rigorous 
statistical-mechanical model to see whether the theoret- 
ically calculated phase diagram includes some of these 
loose-packed structures - the A15 lattice in particular. In 
previous theoretical studies this structure was not con- 
sidered as a trial state of a repulsive colloidal system, 
and the solid part of the phase diagram was [Shared vir- 
tually exclusively by FCC and BCC phasesliErQ We scan 
the phase diagram for the case of square-shoulder soft 
potential, the focus of other studies, and determine the 
range of widths of the soft potential where the A15 lattice 
is the colloidal ground state. 

The paper is organized as follows: in Section II, we de- 
scribe our model of colloidal self-organization and estab- 
lish the analogy between these systems and soap froth. 
We digress to describe the Kelvin problem and its con- 
jectured solution. We show that this analogy leads to 
a second global principle of self-assembly - the primary 
one being the principle of maximum packing entropy - 
and that the two mechanisms give rise to frustration. 
We apply these ideas to a dendrimer compound (which 
crystallizes into the A15 lattice) as well as to other sys- 
tems. In Section III we reexamine the phase diagram 
of the hard-core-square-shoulder interaction, which cap- 
tures many features of real colloidal systems. We analyze 
it within the cellular free volume approximation using a 
numerical model whose main advantage is that the col- 
loidal interaction is treated nonperturbatively. Section 
IV concludes the paper. 



II. COLLOIDAL SYSTEMS AS 
AREA-MINIMIZING STRUCTURES 

Recently, dendrimers composed of a poly(benzid. 
ether) core segments decorated with dodecyl chainsEEl 
were synthesized with the intention of producing a 
molecule with a conical, fan-shaped architectureo Since 
it is known that many dendrimers spontaneously self- 
assemble into supramolecular clusters, the underlying ra- 
tionale was motivated by the possibility of creating spher- 
ical micellar-like objects a few nanometers in diameter. 
These dendrimers do indeed form spheres which, in turn, 
form a crystal lattice with the Pm3n space group, also 
known as the A15, Q 223 , and (3- Tungsten lattice. This 
lattice belongs to the cubic system, and its unit cell in- 
cludes 8 sites which can be divided into 3 pairs of colum- 
nar sites and 2 interstitial sites. The columnar sites lie 
evenly spaced along the bisectors of the faces of the unit 



cell and can be thought of as forming three mutually 
perpendicular and interlocking columns. The interstitial 
sites fill out the space between the columns: one is at 
the center of the cell and the other one is at the vertex 
(Fig- 0). 




FIG. 1. Various lattices: (a) Face-centered cubic, (b) 
body-centered cubic, and (c) A15 lattices. 

Though the structure pf .dendrimer crystals has been 
studied in great detailjl3't3 it remains unclear why 
the dendrimer assemblies should make the A15 lattice. 
The micelles axe nearly spherically symmetric and very 
monodisperse:tj given their chemical composition, the in- 
teraction between the micelles must be predominantly 
steric. If they were of uniform density, the interaction 
arising from the impenetrability of the micelles would 
be well-described by a spherically symmetric hard-core 
potential. In this case, one might expect the spheres to 
assemble into an FCC lattice to maximize their positional 
entropy - because it has the largest packing fraction, it 
maximizes the volume available to each sphere. O-EJ On 
the other hand, the A15 lattice is rather loosely packed, 
with the same packing fraction as that of the simple cu- 
bic (SC) lattice. Thus the A15 lattice is very inefficient 
from the point of view of the center-of-mass entropy of 
each sphere. 

If the micelles were structureless hard spheres, their 
free energy would depend only on their position. How- 
ever, the dendrimers have a well-defined structure con- 
sisting of a more or less compact core of benzyl ether rings 
and a floppy, squishy corona of alkyl chains. In this case, 
the stability of a certain arrangement of micelles does 
not depend only on the hard-core repulsion between the 
cores but also on the interaction between the brush-like 
coronas: the larger the overlap between the neighboring 
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micelles, the more constrained the conformations of the 
chains within the coronas and the smaller their orienta- 
tional entropy (Fig. ||). This effective interaction is, of 
course, repulsive but short-ranged: at distances larger 
than the diameter of the micelles, their coronas do not 
overlap and the interaction vanishes. 



potential 




distance 



FIG. 2. Origin of the soft intermicellar potential: over- 
lap between the micelles reduces orientational entropy of 
chains within the corona (shaded area), which gives rise to 
short-range repulsion. 

Though the hard and the soft part of the repulsion 
both arise from the steric interaction between the den- 
drimers, their dependence on the density is very different. 
The hard cores lead to an inaccessible volume, while the 
matrix of interpenetrating soft coronas leads to a softer 
cntropic repulsion. We can regard the matrix of coronas 
as a bilayer of dodecyl chains wrapped around each hard 
core. The free energy of these bilayers decreases mono- 
tonically with thickness. However, the volume of this 
soft matrix is fixed by the difference between the total 
volume and that of the hard cores and can be written as 
the product of the total area (A) of these bilayers and 
their average thickness (d) , so that at a given density 

Ad = constant. (1) 

Since the soft repulsion of the tails favors larger bilayer 
thicknesses d, the coronal entropy is maximized when the 
area A is minimized, or, in other words, when the inter- 
facial area between the neighboring micelles is a minimal 
surface (Fig. |3|). We should emphasize that the interface 
between the micelles is a mathematical concept that em- 
bodies the membrane-like, two-dimensional character of 
the interdigitated coronas, and it does not correspond to, 
say, the position of a particular segment of the dendrimer 
molecule. 
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FIG. 3. A two-dimensional illustration of the minimal-area 
rule: the hard cores are embedded in the matrix of coronas, 
and the volume of the latter is given by the product of its 
area (dotted line) and the average separation of the cores. 

It is essential to note that this minimal-area princi- 
ple is incompatible with close packing and can favor lat- 
tices other than FCC, thereby giving rise to frustration in 
these systems and a rich phase diagram. If the free energy 
of the lattice of micelles depended only on the interfacial 
energy, the system would behave as periodic membrane 
enclosing bubbles of equal volume - as an ideal, dry soap 
froth. At this juncture the solution to the problem of 
minimizing the interfacial area of a set of equal volume, 
space-filling cells is not known. However, it is known that 
the area of the BCC "foam" is smaller than that of the 
FCC foam and that the A15 foam has a smaller area still. 



A. Kelvin's problem 

The problem of finding the ideal configuration of a 
soap froth was introduced by Kelvin in 1887 while study- 
ing how light propagated in a crystal, the relation being 
based on the now abandoned notion of aether. He real- 
ized that the dry soap froth problem could be cast in the 
mathematically precise form: What regular partition of 
space into cells of equal volume has the smallest area of 
cells? As is almost always the case, the plain and simple 
formulation of the problem is a harbinger of its complex- 
ity. 

Kelvin built on Plateau's investigations of the stability 
of soap films, summarized by two rules which represent 
mechanical equilibrium: In an equilibrium froth, (i) ad- 
jacent faces meet at an angle of 120° and (ii) the edges 
(the so-called Plateau borders) must form a tetrahedral 
angle of 109°28'. It follows that arrangements with more 
than 3 faces meeting at a common edge are unstable as 
are junctions with more than 4 edges. At that time, 
these rules did not have a theoretical background and 
neither did Kelvin's work. Kelvin's approach relied on 
experiments based on soldered wire frames dipped into 
soap solution and the observaticm^of the evolution of the 
soap film spanned by the frame.113 Kelvin was led to the 
conclusion that his problem was solved by a lattice of 
polygons with the topology of an orthic tetrakaidecahe- 
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dron consisting of 6 quadrilateral and 8 hexagonal faces. 
In more modern terms these shapes are known as the 
Wigner-Seitz cell of the BCC lattice (Fig. |). Kelvin 
subsequently worked out the exact shape of the faces and 
showed that to satisfy the Plateau rules, the edges of the 
tetrakaidecahedra must be slightly curved and the hexag- 
onal faces must be somewhat nonplanar. 

Though a conjecture, Kelvin's partition of space was 
thought to be the solution of the problem by the math- 
ematics community, despite the realization the proof 
might be highly elusive.EII Indeed, these types of prob- 
lems are notoriously difficult. For example, even the 
Plateau rules themselves were nothing but experimen- 
tal facts until 1976 when_they were put on a firm theo- 
retical footing by TayloiO and it was only in 1999 that 
Hales proved that the regular hexagon is the solution of 
the two-dimensional variant of the Kelvin problem.CJ We 
also note that the related Kepler problem of packing hard 
spheres as densely as possible turned out to be equally 
challenging: the well-known FCC or HCP packing was 
demonstrated to be the most efficient in 1998, the proof 
being furnished again by Hales £3 
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FIG. 4. Rhombic dodecahedron/FCC lattice (a), Kelvin's 
tetrakaidecahedron/BCC lattice (b), and Weaire-Phelan min- 
imal surface/A15 lattice (c). 



Kelvin's conjecture stood unchallenged for more than 
a century until 1994 when Weaire and Phelan discovered 
that a froth with the symmetry of the AL5, lattice has 
an area smaller than BCC lattice by 0.3%.Ea While this 
may seem to be a small difference, it is, in fact, signifi- 
cant: the relative difference in the areas of Kelvin's BCC 
structure and the FCC-type partition of rhombic dodec- 
ahedra (unstable as a froth because it contains vertices 
joining 8 edges) is 0.7%. Thus the A15 foam is a 50% 
improvement on this scale. 

The A15 foam consists of 6 Goldberg tetrakaideca- 
hedra, each with 2 hexagonal and 12 pentagonal faces, 
which form three sets of interlocking columns, and 2 ir- 
regular pentagonal dodecahedra at the interstices. Since 
it is composed of two types of cells, the Weaire-Phelan 
soap froth is obviously less symmetric than Kelvin's 
froth. However, it was not surprising to some that it 
has a rather small surface area. The A15 lattice is de- 
rived from the polytope {3, 3,5}, a partition of positively 
curved space consisting of 120 regular dodecahedral bub- 
bles. As suggested by Kleman and Sadoc, positive spatial 
curvature relieves the frustration brought about by the 
incompatibility of the optimal local arrangement of bub- 
bles and the structure of flat three-dimensional space. 
The polytope {3,3,5} should be regarded as the ideal 
template whence the various tetrahedrally close-packed 
(TCP) lattices (the layered TPC lattices are also known 
as Frank-Kasper phases) are derived via decurving, i.e., 
by substituting some of the dodecahedral cells by bub- 
bles with 14 (= 12 pentagonal + 2 hexagonal), 15 (= 12 
pentagonal + 3 hepamraal) , or 16 (= 12 pentagonal + 
4 hexagonal) faces. EjESEj Cells with 13 faces are forbid- 
den for topological reasons, and-those with more than 16 
faces are dynamically unstableO 

There are 24 known ways of decurving, the polytope 
{3, 3, 5} and thus 24 TCP crystal latticesB With 2 types 
of bubbles and 8 bubbles per unit cell, the A15 lattice is 
among the simplest: the unit cell of the most complicated 
structure, the so-called I lattice, consists of 228 bubbles 
that include all 4 types of bubbles &3 Given the rationale 
for the success of the A15 lattice, it is not unreasonable to 
expect that other TPC lattices may have an even smaller 
surface area. Indeed, the discovery of Weaire and Phe- 
lan renewed the interest in the field. Subsequently some 
of the remaining TCP-type soap froths have been stud- 
ied,^ facilitated by Surface Evolves- a remarkable soft- 
ware package developed by Brakke£3 In addition, other 
classes of periodic and quasiperiodic partitions based on 
Kelvin's BCC and Williams' body.=centered-tetragonal 
(BCT) bubbles have been suggested.Ej At this time, how- 
ever, the A15 foam has not been bested - it stands as the 
tentative solution of Kelvin's problem. 

The connection between the observed crystal structure 
and the soft part of the intermicellar potential should 
now be clear: if it were absent, the hard cores would 
favor an FCC lattice, but instead the equilibrium struc- 
ture is a different lattice with a larger bulk free energy 
but a smaller surface free energy. The only variables in 
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this model are the parameters of the surface interaction, 
which we will estimate roughly. In view of the experimen- 
tal data on soft spheres, we will compare the FCC, BCC, 
and A15 lattices and show that in the case of the den- 
drimer micelles, it is reasonable to expect that the area- 
minimizing A15 lattice is favored over the close-packed 
FCC lattice. 



B. Bulk free energy 

Within the framework of the two competing ordering 
principles, we propose a simple and approximate theory 
of colloidal crystals where the bulk and the surface terms 
are coupled only through the constraint of fixed volume. 
As far as the bulk free energy is concerned, the micelles 
will be treated as hard spheres of diameter a, and the 
surface interaction will be calculated as if the matrix of 
coronas were a thin structureless layer so that we may 
neglect the effects of curvature. 

We start our analysis with the bulk term. The config- 
uration integral of a one-component classical system of 
N particles confined to a volume V is 

i r N 

Z =^Wm / II dr * exp(-C/(r 1 ,r 2 ,...r Ar )/fc B T), 



(2) 



where A = yjh 2 /2-Kmk B T is the thermal de Broglie wave- 
length and {/(r 1) r 2 , . . . r N ) = \ ufarj) is the to- 
tal interaction energy consisting of pairwise interactions 
u(rj,rj) = uflrj— r,-|). Despite the simplicity of the hard- 
core potential, Z cannot be calculated analytically, and 
one must resort either to numerical approaches, such as 
Monte Carlo simulations, or to approximate analytical 
methods. The latter are more appropriate for our pur- 
poses since we are seeking a simple, heuristic explanation 
of crystal structure. When considering crystal phases we 
can assume that the particles are localized within cells 
formed by their neighbors and thus the configuration in- 
tegral breaks up into N single-particle integrals so that 
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where Vo is the volume of the cell (usually the Wigner- 
Seitz cell), u c g(r) is the effective potential felt by the 
particle, and N\ has been absorbed by factorization of 
the partition function. This approximation is valid only 
for lattices where all sites are equivalent; if the unit cell 
of the crystal consists of inequivalent sites, Eq. (j^) can 
be embellished accordingly. 

As natural as this approximation may seem for a pe- 
riodic arrangement of particles in a crystal, it neglects 
any correlated mo tion .of neighbors and the associated 
communal entropyOCj Nevertheless, the cellular model 



often gives quantitatively accurate results that are in 
good agreement with more complete numerical simula- 
tions. For hard spheres, where the cellular theory be- 
comes exact in the high-density limit, t he ag reement with 
Monte Carlo results is in fact excellent In this case, 
each particle is assumed to be uniformly smeared over 
its reduced Wigner-Seitz so that it cannot overlap with 
its neighbors. As a result, the effective potential is quite 
simple: within the volume the center of mass is allowed 
to trace out (the so-called free volume) and infinite oth- 
erwise. Thus 
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where Vf is the free volume whose shape reflects the sym- 
metry of the lattice.Ej We note that temperature drops 
out of this problem and that the interaction is entirely 
entropic. 

Within the cellular free- volume theory, the free energy 
of hard spheres becomes a purely geometrical issue. All 
one has to do is to calculate the volume of the Wigner- 
Seitz cell after a layer of thickness a/2 has been peeled off 
of its faces, which is the volume accessible to the particle's 
center of mass. In the FCC lattice, the free volume has 
the shape .of a rhombic dodecahedron just as the Wigner- 
Seitz cell.E 2 ] On the other hand, in the BCC lattice it 
remains an orthic tetrakaidecahedron only at rather low 
densities far below the freezing point. At higher densities, 
where the hard spheres form a solid phase, the square 
faces become absent rendering the free volume a regular 
octahedron. For the FCC and BCC lattices, the bulk free 
energy per dendrimer micelle reads 



F^ ulk = -k B T In 
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where X stands for FCC or BCC, n — pa 3 is the reduced 
number density (a being the hard-core diameter), and the 
term — 3k B Thi(a/ X) which is independent of the lattice 
structure has been dropped. The coefficients a FCC = 
2 5 / 2 and a BCC = 2 2 3 1/2 reflect the shape of the free 
volume, and (3 FCC = 2 1 / 6 and /3 BCC = 2- 2 / 3 3 1 / 2 specify 
their size. 

The free volume of the A15 lattice, which includes 
two types of sites, is a bit more complicated. As deter- 
mined by the Voronoi construction subject to the con- 
straint that all cells have equal volume, the shapes of the 
free volumes are irregular pentagonal dodecahedra and 
tetrakaidecahedra with two hexagonal and twelve pentag- 
onal faces rather than regular polyhedra. Their volumes 
cannot be expressed in an amenable analytical form, and 
thus we have calculated them numerically for a range of 
densities. However, we note that the exact free volumes 
can be approximated very well by replacing the dodec- 
ahedra and tetrakaidecahedra by spheres and cylinders, 
respectively (Fig. ||). 
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FIG. 5. Approximate analytical model of the bulk free en- 
ergy for the A15 lattice: free volumes of columnar and inter- 
stitial sites are replaced by cylinders and spheres, respectively. 

This ansatz takes into account that the free volume of 
a columnar site changes anisotropically with density and 
approaches a flattened shape in the close-packing limit, 
whereas the shape of the free volume of interstitial sites 
does not depend on density. We also introduce two ad- 
justable parameters that quantify the fact that the actual 
volumes of the Wigner-Seitz cells are larger than the vol- 
umes of spheres and cylinders, which leave empty voids 
between them. Given that the ratio of columnar and in- 
terstitial sites is 3:1, this leads to the average bulk free 
energy per micelle of 




This formula best agrees with the numerical results for 
S = 1.638 and C — 1.381, where the relative deviation 
from the true bulk free energy is below 0.1% at densities 
higher than n ss 0.8 and not significantly larger at lowei; 
densities. Fg$g, F^, and are plotted in Fig. §B 
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FIG. 6. Bulk free energies of hard-core particles arranged in 
FCC, BCC, and A15 lattice as calculated with the free- volume 
theory. Solid lines correspond to Eqs. (|^) and (ffl), and circles 
are numerical results. 




FIG. 7. Structure of dendrimers that form the A15 lat- 
tice: 1st generation generation dendrimer consists of a core 
segment and three dodecyl chains, 2 nd generation dendrimer 
includes three 1 st generation dendrimers attached to a core 
segment, and so on. The length of the core segment and the 
dodecyl tail is / corc = 0.6 nm and /corona = 1-4 nm, respec- 
tively. Shown here is a 2 nd generation dendrimer. 

C. Surface free energy 

Having calculated the bulk free energy for the three 
lattices, we now turn to the surface free energy, which 
requires a specific model of the soft interaction between 
the particles. To determine which model is the most ap- 
propriate, we should first find out whether the overlap of 
the neighboring particles is large or small. At this point, 
the analysis becomes somewhat less general, because the 
overlap differs from one system to another. Being in- 
terested primarily by the stability of the loosely-packed 
crystal lattices, we new focus on the dendrimers that 
form the A15 crystal.li-3 The relevant quantitative data 
include measurements of the radii of dendrimer molecules 
and the lattice constants of the micellar crystals. Since 
the dendrimers consist of 2, 3, or 4 generations of branch- 
ing benzyloxy segments crowned by dodecyl chains with 
bare radii of <7barc/2 = 2.6, 3.2, and 3.8 nm, respectively, 
we can deduce that the lengths of the benzyloxy core 
segment and the dodecyl chain are i core = 0.6 nm and 
Corona = 1.4 nm (Fig. [?]). The effective diameter of the 
micelles cr e ff can be calculated from the lattice constant 
a. According to Fig. [I], a e g = a/2. This gives a e g = 3.4, 
4.0, and 4.2 nm. Obviously, cr c ff is considerably smaller 
than (Tbare for all generations. Actually, the most con- 
clusive information can be extracted from 4 th generation 
data: 2 nd and 3 rd generation micelles probably have an 
empty center and their true diameter is most likely larger 
than (Tbare- These data show that the effective diam- 
eter of the hemispheric 4 th generation micelles, 4.2 nm, 
does not exceed the diameter of the their benzyloxy core, 
8^corc = 4.8 nm. Thus we can conclude that the hard-core 
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diameter of the monodendrons must be smaller than the 
diameter of the benzyloxy core. Thus not only is there 
a considerable overlap between the coronas, but also the 
dodecyl chains penetrate into the core itself. This indi- 
cates that the interfacial effects, related to the limited 
orientational entropy of the chains, are important in the 
dendrimer system. 

In the absence of a quantitative insight into the inter- 
micellar potential - such as a direct measurement of the 
interaction via optical trappingcij - we model the interac- 
tion of the interpenetrating dodecyl chaijas. as the inter- 
action between grafted polymer brushes." In the high- 
interdigitation limit, which is certainly applicable in our 
case, the free energy of a compressed brush reduces to the 
excluded- volume repulsion of the chains. An argument 
in the spirit of Flory theory gives 
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where £ is a parameter with the dimension of length, 
which determines the strength of repulsion, Nq is the 
number of chains per micelle, and h, the thickness of 
the single corona, is half the average thickness of the 
interdigitated matrix of the chains d. This approxima- 
tion neglects some details of the actual interaction, most 
notably the curvature of the dodecyl brushes in the coro- 
nas. However, these effects are subdominant when the 
dendrimers are packed very closely. In fact, it is known 
that the density is,-«early constant in the volume occu- 
pied by the chains, til and we thus expect that our model 
should provide a robust description of the system. 

The thickness d of the coronal matrix depends on the 
density since the dodecyl bilayer must fill the space be- 
tween the hard cores. lij Thus 
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where Am is the interfacial area per micelle. Since the 
area is proportional to the square of the lattice constant, 
the surface free energy per micelle reads 
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where ^ x is the coefficient determined by the symmetry 
of the lattice and defined by Am = j x a 2 n~ 2 / 3 . Lower 
values of j x correspond to more area-efficient partitions. 
The typical magnitude of this coefficient is set by the 
simple cubic (SC) lattice, which gives 7 SC = 6. The 
more efficient area-minimizing lattices have smaller val- 
ues of 7 X : 7 F cc = 2 5 / 6 3 = 5.345, 7 BCC = 5.306, and 
7 A15 = 5.288.E3 7 BCC and 7 A15 can onk be computed 
numerically, e.g., using Surface Evolver.EZl As a compar- 
ison, the ultimate lower bound for 7 corresponds to a 
sphere (which, of course, is not a space-filling body): 

7 sphcrc _ 22/332/3^1/3 _ 4 836. 

Combining Eqs. (|j), ([?]), and (^), we arrive at a single- 
parameter (£) free energy of micellar crystals: 



Before we calculate the minimal strength of the soft re- 
pulsion necessary to stabilize the area-minimizing struc- 
tures, we need to determine the density of the crystal. 
Not knowing the hard-core radius of the micelles, we 
can only provide an order-of-magnitude estimate of £, 
and for this purpose, it is sufficient to note that unless 
the observed A15 lattice approaches the packing limit 
of n = 1.0, its bulk free energy per micelle is roughly 
lk B T larger than and 2k B T larger than F^ . 

To be concrete, we compare the total free energies at 
n = 0.95, which is well within the high-density regime 
but not quite at the close-packing limit and thus consis- 
tent with the structural data. At this density, the BCC 
lattice becomes more favorable than FCC at I > 0.05cr 
whereas the BCC-A15 transition occurs at I w 0.15er. 
Given that there are iVo = 162 dodecyl chains per 3rd 
and 4th generation micelle, these values of I correspond 
to an entropy of about 0.5fc B and 1.5fc B per chain, respec- 
tively. In other words, this means that if the overlap of 
the micelles is so large that the decrease of orientational 
entropy of a chain due to interdigitation reaches 0.5fc B 
and 1.5/c B , the differences between the surface entropies 
of the BCC and FCC lattices and the A15 and BCC lat- 
tices overcome the corresponding differences of the bulk 
entropies and thus favor the area-minimizing structures. 

These values are physically reasonable and of the cor- 
rect order of magnitude. When unrestricted by other 
chains, each chain has a few orientational and confor- 
mational degrees of freedom. Thus its entropy is a few 
k B , and so it can easily loose 0.5fc B , 1.5fc B , or more en- 
tropy upon interdigitation. If we plug these numbers 
back into the free energy, we learn one more thing about 
these systems: their energetics is controlled primarily by 
the surface term, which is a direct consequence of the 
three center-of-mass degrees of freedom associated with 
the hard cores being greatly outnumbered by the several 
hundreds of internal degrees of freedom associated with 
the dodecyl chains. In other words, in the dendrimer 
system the "squishiness" of the particles wins over their 
hard cores. 



D. A new paradigm 

Our proposed model provides a novel way of looking 
at the self-organization of colloidal crystals. By comple- 
menting the well-known close-packing rule, which con- 
trols the stability of hard spheres, with the minimal-area 
rule, which stems from the additional short-range repul- 
sion between the particles, we have shown that the equi- 
librium ordered structures are a result of frustration be- 
tween two incompatible requirements. This picture im- 
plies that if the corona is thin compared to the core, the 
colloids will behave as hard spheres and form a close- 
packed lattice such as FCC, but as it grows thicker, an 
area-minimizing structure - the A15 lattice - should be 
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observed. In between the two extremes, there may be a 
spectrum of lattices that neither maximize the packing 
fraction nor minimize the interfacial area but represent 
a reasonable compromise for a given intermicellar poten- 
tial. However, in some systems coexistence of the FCC 
and A15 lattices could be energetically preferable to an 
intermediate structure such as BCC, provided that the 
density is not too high to destabilize the A15 lattice. 

This observation is consistent with the experimentally 
determined structures found in other colloidal materi- 
als, such as crystals of aurothiol particles consisting of a 
gold crystallite core and covered by about 50 n-alkylthiols 
where n = 4, 6, or 12. This system is remarkably close 
to the dendrimer system in both size and structure. The 
diameter of the metallic cores can be varied from 1.6 to 
3.1 nm, and the length of the (fully extended) coronal 
chains is between 0.6 to 1.56 nm. The alkylthiol chains 
were adsorbed to the gold core with the sulfur atom, leav- 
ing the outer part of the corona chemically identical to 
the dendrimer coronas. Depending on the relative size 
of the corona with respect to the diameter of the core, 
either the FCC or BCC lattice was observed. In addi- 
tion, in some samples the BCT lattice was found as a 
moderately anisotropic variation of the BCC lattice with 
c/a ~ 1.15. (We will see later on that this lattice is not 
really unexpected although it departs somewhat from our 
model: in a BCT crystal, both the bulk and surface en- 
tropy are larger than their BCC counterparts.) 

A particularly interesting feature of this system is that 
the size of the core can be varied continuously, so that 
the FCC-BCC transition can be located very precisely. 
For the range of core radii explored in the study,L3 an 
FCC-BCC transition was found in particles covered with 
hexylthiol chains: the FCC lattice is stable if the ratio of 
the thickness of corona and the core radius is smaller than 
about 0.73, and BCC or BCT lattices are observed oth- 
erwise. The structural parameters of the aurothiol parti- 
cles were measured in some detail and if we identify the 
hard core with the gold nanocrystal, we find that the re- 
duced density at the transition is about jt-fcc-bcc = 0.4. 
Proceeding along the same lines as before, this gives an 
entropy decrease of about 0.05fc B for each of the approx- 
imately 150 chains in the corona. Although still reason- 
able, this estimate is by an order of magnitude smaller 
than in the dendrimer case. However, we note that be- 
cause of the high coverage of the gold cores with the ad- 
sorbed alkylthiols, the effective hard-core diameter of the 
particles is most likely larger than the gold core diameter. 
In this case, the reduced density at the transition would 
be larger, and this would lead to a larger entropy decrease 
per alkylthiol chain. Nonetheless we can still combine 
this figure with its dendrimer counterpart to bound the 
value of the entropy per coronal chain at the FCC-BCC 
transition between 0.05 and 0.5fc B , which can serve as an 
estimate for future studies. 

Similar behavior was discovered in polystyrene- 
polyisoprene diblock copolymers dispersed, ah decane, a 
solvent preferential for the polyisoprene £313 These di- 



blocks spontaneously form micelles with a polystyrene 
core and polyisoprene corona. Indeed, micelles based on 
diblocks with a core segment containing from 1.5 to 2 
times as many monomeijs-as the coronal segment crystal- 
lize into an FCC latticeO On the other hand, in copoly- 
mers made of blocks with the same number of monomers 
the BCC lattice was observed. Although the A15 lattice 
was not seen in this system, we conjecture that it could 
be found in copolymers with the polyisoprene block suf- 
ficiently longer than the polystyrene block. 

While rare in crystalline systems, we mention that 
the A15 structure is not uncommon in lyotropic sys- 
tems, typically containing lipids dispersed in water ma- 
trix. These systems self-arrange in either direct or in- 
verted micelles, -which are known to form a variety of 
cubic structures.E3cil However, the existence of the A15 
lattice in lyotropic systems should not be as surprising 
as in colloidal systems. The cohesive force of these mi- 
celles is the hydrophobic interaction which means that 
their size is not as well-defined as in colloids. If there can 
be several types of micelles the A15 lattice is not hard 
to assemble. The other difference with respect to col- 
loids is that the water matrix that encloses the micelles 
is truly fluid and thus much more similar to actual wet 
soap froths with relatively large liquid content. At the 
same time, in the lyotropic systems the effective inter- 
micellar potential may not be necessarily dominated by 
steric effects, implying that the stability of a particular 
ordered micellar structure can depend on the chemical 
composition of lipids. 

Specific effects like these can not easily be incorpo- 
rated into our framework: our goal was not to develop 
an elaborate description that would cover the many de- 
tails that are captured in more or less involved theories, 
such as self-coBsktent field theoryH and various Monte 
Carlo schemes .I13e3 Instead, we have proposed a heuris- 
tic model that describes the essential physics of many 
colloids and clearly exposes the frustration that is in- 
troduced by dressing hard spheres with a soft repulsive 
interaction. This model provides an explanation of the 
results of the more rigorous theories with the added ad- 
vantage that it is very simple, yet can give reasonable 
semiquantitative predictions about the stability of the 
different lattices. Obviously, it can be improved by intro- 
ducing a more refined model of the interfacial interaction 
which would account for the curvature of the coronal ma- 
trix, the strain of the coronas into the interstitial regions, 
solvent effects, and related phenomena. 

As useful as it might be, like all coarse-grained ap- 
proaches, the proposed theorjil3 has some limitations. 
Our model is tailored for a broad yet specific class of col- 
loidal particles that interact by short-range potentials, 
and it may not be suitable for systems characterized by 
long-range forces, such as unscreened or partly screened 
electrostatic forces. The range of potential is determined 
by the requirement that the nearest-neighbor interaction 
be much stronger than the interaction with the rest of 
the particles; otherwise we could not define the concept 
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of the interface. Additionally, we have decoupled the 
bulk and interfacial free energies, or the hard and soft 
parts of the potential. We suspect that a limit of the full 
problem exists in which the minimum-area principle is a 
mean-field approximation. 

In order to gain insight into colloidal crystal struc- 
ture and to strengthen the case for our new principle, we 
will, in the next Section, calculate the phase diagram for 
a certain interparticle potential, and identify the range 
of stability of the A15 lattice using a "single-particle" 
Monte Carlo algorithm. The significance of these results 
is twofold: at the qualitative level, they suggest that the 
phase diagram of soft spheres can include a variety of 
solid phases, not just FCC and BCC lattices, and at the 
quantitative level they will provide guidance for a more 
complete numerical analysis of the system. 

III. SHORT-RANGE POTENTIALS AND 
SOLID-SOLID PHASE TRANSITIONS 

A. Crystal architecture 

The FCC lattice, which, as we know now for certain,El 
is the closest regular packing of hard spheres with cu- 
bic symmetry, with each site enclosed completely by the 
12 equidistant neighbors. In the case of a purely hard 
core interaction, the free energy consists solely of the en- 
tropic contribution and is a convex function of volume. 
To model the effect of the soft coronas in the systems we 
are discussing, we add to the hard core a soft potential 
of range 5. Note that we are replacing an entropic in- 
teraction with an energetic term which should make no 
difference to our argument but is more convenient in our 
calculational scheme. At densities low enough so that 
the average distance between the particles is much larger 
than 6, the overlap of the coronas is small and so the 
free energy is dominated by the entropy, as if the extra 
potential were not there. On the other hand, in the high- 
density limit where the separation between the hard cores 
is smaller than 6, each particle feels the soft potential of 
all neighbors. If the potential is flat enough, the average 
field may in fact not change dramatically with position 
and then the particle will behave essentially as if it were 
a hard sphere moving in a constant potential, its free en- 
ergy being given by the entropic contribution of the hard- 
core interaction shifted by the energy of the overlapping 
coronas. For certain forms of the soft potential, the tran- 
sition from the low-density to high-density regime can be 
rather abrupt (Fig. |J). Upon compression, the free en- 
ergy of the colloid will then change quite sharply and 
may no longer be a convex function of volume. In this 
case, the low-density and high-density behavior will be 
separated by a region of coexistence between expanded 
and condensed FCC lattices, implying an isostructural 
transition between thcm.cl 



Coexistence between the expanded and condensed 
FCC phases is not the only possibility: other crystal lat- 
tices may exist in this density regime where it may be 
favorable for the spheres to be configured in a structure 
with a somewhat higher entropy but lower energy. For 
example, the BCC lattice has a more open architecture 
compared to FCC where the "cage" of each particle is 
defined by the 12 nearest neighbors. In the BCC lat- 
tice, the 8 nearest neighbors do not enclose the particle 
completely: as the Wigner-Seitz cell shows us, even in 
the hard-core limit, each particle also interacts with the 
6 next-to-nearest neighbors. At any fixed density, the 8 
neighbors from the first coordination shell are closer to 
the particle than in the FCC crystal, but the 6 second- 
shell neighbors are further away. As a result, BCC hard 
spheres have a lower entropy than the FCC hard spheres, 
but the additional soft potential can stabilize the BCC 
lattice. The more structured arrangement of sites implies 
that upon compression, the energy of the BCC crystal 
increases more gradually than its FCC counterpart: the 
central particle first overlaps with the first-shell neigh- 
bors alone and in this density regime the BCC free en- 
ergy can be lower than the free energy of the coexisting 
expanded and condensed FCC phases, provided that the 
soft potential overlaps with all 12 neighbors of the FCC 
lattice. At somewhat higher densities the particles in 
the BCC lattice interact with all of the 14 neighbors and 
then both the entropic and energetic terms disfavor the 
BCC lattice. Similarly, the stability of the A15 lattice 
and other more structured phases is facilitated by their 
even smoother, less step-like free energy as a function of 
volume. 

This view of the energetics is closely related to the 
structure of the pair distribution function of a particular 
lattice. As long as the interaction between the particles 
is short-range, one only has to consider the pair distri- 
bution function at short separations, which directly re- 
flects the shape of the corresponding Wigner-Seitz cell. 
If the distribution function has a single peak as in the 
highly symmetric FCC crystal, a soft short-range repul- 
sion between the particles can destabilize the lattice with 
respect to other lattices with more complex pair distribu- 
tion functions. As noted above, these include the BCC 
lattice with two peaks at short distances (Fig. ^ but 
also the anisotropic variants such as the BCT latticet£l 
where the 6 BCC second-shell sites are split into a sub- 
sheU-pf 4 sites and a subshell of 2 sites, and the BCO lat- 
ticetJ where the fourfold BCT next-to-nearest-neighbor 
subshell is further split in two subshells. (Curiously, the 
anisotropic version of the FCC lattice - the face-centered 
orthorombic (FCO) structure has not yet been observed 
in a colloidal or related system.) 
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FIG. 8. Schematic of the free energy profile for a model 
hard-core-soft-shoulder potential. The short range of the 
shoulder makes the free energy a nonconvex function of vol- 
ume, and thus induces either isostructural or structural tran- 
sitions in the solid phase. Solid line: FCC; dashed line: BCC; 
dotted line: A15. 

This argument can be extended to include other, less 
symmetric lattices. An example is the A15 lattice. By 
considering the pair distribution function, we see that, on 
average, each site interacts with 13.5 neighbors: through 
the faces of their Wigner-Seitz cells, the columnar and 
interstitial sites interact with 14 and 12 neighbors, re- 
spectively, and there are 3 times as many columnar sites 
as the interstitial sites. Of these 13.5 neighbors, 1.5 par- 
ticles are closest to the central particle, followed by two 
shells of 6 particles at larger distances (Fig. [)]). While it is 
obvious that this lattice must have a smaller entropy than 
BCC, for a suitably chosen soft potential it could have 
a lower net free energy. An extreme case of the tradeoff 
between entropy and energy is the diamond phase ob- 
served in a numerical study of packing of star poijfmers 
which interact with an extremely soft potential.113 The 
diamond lattice with 4 nearest and 12 next-to-nearest 
neighbors is clearly unfavorable from the entropic point 
of view, but apparently the pair distribution function is 
so strongly peaked at next-to-nearest neighbors that for 
certain types of short-range interparticle potentials the 
energy of the diamond lattice is very low. 

While these ideas provide additional perspective and 
complement our foam model described above, it is not 
completely clear whether the structural characteristics 
captured by the pair distribution function can be easily 
related to the minimal-area principle and Kelvin prob- 
lem. We will not pursue this interesting question here; 
instead, we now turn to statistical mechanics of soft 
spheres. 

B. Phase diagrams 

The phase behavior of soft-sphere systems has had on- 
going study over the years. A variety of purely repul- 
sive forces have been considered, including power law 
potentials as the high-temperature limit of the van der 



Waals interaction,!^ screened Coulomb potentials withcJ 
or without a hard core,™ and square shoulder poten- 
tials&L}E3 as well as their variants. c3 With its flat plateau 
and as sharp a cutoff as possible, the square-shoulder 
potential is very sensitive to the structure of the pair 
distribution function and thus the quintessential short- 
ranged interaction. This is the reason why it is used so 
often to analyze solid-solid transitions. We shall use the 
square-shoulder interaction in the following, in order to 
provide a comparison with earlier studies. 
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FIG. 9. Structure of the FCC (solid line), BCC (dashed 
line), and A15 lattice (dotted line) represented by the radii 
of a few nearest coordination shells weighted by the corre- 
sponding coordination number M. (If peaks were multiplied 
by (5-functions, this would be the T — pair distribution 
function integrated over the solid angle.) FCC is the most 
closed lattice, each particle having 12 equivalent neighbors, 
whereas BCC and A15 are more open structures with 2 and 
3 inequivalent neighbor sites, respectively. 

The hard-core square-shoulder (HCSS) system is char- 
acterized by a pairwise potential of the form 

!oo r < (T 
e a < r < a + S (11) 
r > 0- + S 

where S and e > are the thickness and the height of the 
shoulder, respectively (Fig. |l^). This system has already 
been studiecfrca by using density-functional perturbation 
theory, which treats the behavior of the system primarily 
by the hard-core interaction and treats the additional 
soft potential as a perturbation. If the free energy and 
the structure of the reference state (F re {) are known, the 
total free energy rpay be expanded to linear order in the 
perturbative termE2: 

F[p(r)} = F xei [p(r)} + ~ j dr x dr 2 ^(iq, r 2 ) 0(r ia ), 

(12) 
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where r = T2 — T\ describes the relative position of the 

(2) 

two particles, p rc{ (ri, is the pair distribution function 
of the reference state, and pO rt is the perturbative part 
of the potential. The inputs of this scheme are the pair 
distribution function and the free energy of the reference 
system. Typically, the free energy of the solid HC phases 
are described within the so-called cellular free volume 
theoryoO and the corresponding pair distribution is cal- 
culated from the Gaussian one-particle density, whereas 
the fluid free energy and pair distribution are described 
by the Carnahan-Staciing formula and the Verlet-Weis 
formula, respectivelyJlj 

Within perturbation theory, it has been found that the 
solid part of the phase diagram can indeed be very com- 
plex, and that its topology depends very delicately on 
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FIG. 10. Hard-core-soft-shoulder potential 

the parameters of the potential. For relatively narrow 
shoulders, the phase diagram differs from its hard-core 
counterpart only in the expanded FCC-condensed FCC 
transition that occurs at densities near the close-packing 
density and terminates at a critical point. As the shoul- 
der becomes broader, the region of coexistence between 
the expanded and condensed phases shifts towards lower 
densities. Eventually, the expanded FCC structure is re- 
placed in part by the BCC at intermediate temperatures, 
i.e., at temperatures high enough so that the shoulder 
does not appear very high, yet low enough so that it is 
not irrelevant. According to this analysis, the transition 
between the fluid and condensed FCC phase can be ei- 
ther direct or indirect, the intervening phases being BCC, 
expanded FCC, or both. Similar results were obtained 
for a sloped shoulder and a square shoulder followed by 
a linear ramp, except that these potentials give rise to 
several new features of the phase diagram, such as the 
expanded BCC-condensed BCC transition and a triple 
point. These predictions are consistent with the limited 
Monte Carlo results available,B but most of the features 
remain to be verified by rigorous numerical studies. 

The perturbation analysis has provided valuable infor- 
mation on the behavior of HCSS spheres, but its pre- 



dictive power is limited to high temperatures where the 
expansion is valid. Although some additional informa- 
tion can be obtained by interpolating between the high- 
temperature and zero-temperature data, this leaves out 
the most interesting part of the phase diagram. A natural 
extension of this approach is to calculate the configura- 
tion integral nonperturbatively. This can be done, for ex- 
ample, within a mean-field approximation in the cellular 
model where each particle is assumed to move indepen- 
dently in the average field of its neighbors. Such an ap- 
proximation can be used within the cellular model by as- 
suming an appropriate ansatz for the probability density 
- usually a sum of Gaussians - which then leads to a self- 
consistency relation determined by the requirement that 
probability densities of all particles be the same, thereby 
fixing the parameters of the ansatzn Here we follow the 
spirit of this approach but we note that in the case of 
a system interacting with flat potential such as HCSS, 
the Gaussian ansatz may not be adequate. Moreover, in 
the A15 lattice the probability distribution of particles in 
columnar positions should be very anisotropic and in fact 
nonspherical. Given that our main task here is to exam- 
ine the possibility that the phase diagram of the HCSS 
system includes the A15 lattice, it seems appropriate to 
capture the geometrical details of the different structures 
by calculating the configuration integral numerically. 

Our procedure places the central particle within a cage 
of its neighbors that are themselves allowed to sample a 
distribution of positions such that the probability densi- 
ties of all particles at equivalent sites are the same. This 
is achieved by sampling the actual probability density of 
the central particle continuously, and recreating identical 
distributions at the equivalent neighbor positions.E3 Af- 
ter equilibrating the system for a certain number of steps 
and checking for the consistency of the distributions of 
the central particle and the neighbors, we calculate the 
configuration integral and the free energy. This approach 
could be dubbed a "single-particle Monte Carlo method" ; 
we note, however, that it is in fact a numerical variant of 
the self-consistent mean-field approximation and a rela=. 
tive of the variational method proposed by KirkwoooEJ 
rather than a simplified variant of a true multiparticle 
Monte Carlo integration. Our procedure is single parti- 
cle because we approximate the multiparticle probability 
density by a product of single-particle probability densi- 
ties and thus ignore long-wavelength excitations. 

To ensure that the model reduces to the quantita- 
tively successful cellular free-volume theory used in Sec- 
tion Hell, we modify the HCSS interparticle potential by 
not allowing the hard core of the particle to leave the cell. 
This gives the correct free energy for small shoulders, 
(e "nrkpT) in good agreement with Monte Carlo analy- 
sesflH as well as for large shoulders, where the particles 
behave essentially as hard spheres moving in a constant 
potential. We expect that between these two limiting 
cases, the correlations between nearest neighbors (and 
thus the free energy) should also be well approximated. 
While this method can not include the fluid phase, we 
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are really interested in the stability of the A15 lattice so 
that the only relevant aspect of the fluid-solid transition 
is melting. Instead of going into the details of the fluid 
HCSS phase and analyzing the transition with, for ex- 
ample, Ramakrishnan- Yussouff theoryta we can estimate 
the melting curve by extrapolating the zero-temperature 
FCC melting point, which is, as we will see, far enough 
from the other phase boundaries to suggest that the lat- 
ter persist within the temperature range that we study. 
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FIG. 11. Phase diagram of the hard-sphere-soft-shoulder 
potential for various shoulder widths: (a) 5 = 0.2 and (b) 
8 = 0.25. For these widths there is no expanded BCC lattice. 
The A15 lattice appears as one of the intermediate struc- 
tures, interposed between the expanded and condensed FCC 
structures. The A15 islands of stability are elongated along 
the temperature axis, indicating that this structure should 
be sought after by varying the density rather than temper- 
ature. Note that the single-phase regions are small com- 
pared to coexistence regions. The dashed lines indicate the 
zero-temperature expanded FCC-fluid (F) coexistence. 

Using this approach, we have analyzed the solid part 
of the phase diagram at low temperatures, focusing on 
the FCC, BCC, and A15 lattices. The coexistence be- 



tween the phases was determined using the Maxwell 
double-tangent construction, and the results are shown in 
Figs, [ll] and |l| for shoulder widths S/a = 0.2,0.25,0.3, 
and 0.35. 



k B T/e 




1.2 1.4 



k B Tle 




FIG. 12. Phase diagram of the hard-sphere-soft-shoulder 
potential for various shoulder widths: (a) 8 — 0.3 and (b) 
8 = 0.35. For these widths there is no condensed BCC lattice. 
As in the previous figure, the A15 lattice appears as one of 
the intermediate structures, interposed between the expanded 
and condensed FCC structures. 

The main features of the phase diagram can be sum- 
marized as follows: 

(i) We find that the A15 lattice can be a stable state of 
the HCSS system between the expanded and condensed 
FCC phase. As expected from the structure of the pair 
distribution function, this lattice is stable for shoulders 
neither too narrow (which would make the particles too 
similar to ordinary hard spheres) nor too broad (which 
would destroy the comparative advantage of the A15 
lattice over the FCC and BCC lattices). The minimal 
and the maximal shoulder widths roughly correspond to 
8 = 0.2 and 0.35, respectively. Note that these widths 
are consistent with the structural parameters of the den- 
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drimer compound where this lattice was observed exper- 
imentally, although the effective hard core of the micelles 
may not necessarily coincide with the benzyl inner seg- 
ment of the dendrimers. 

Temperature is also a crucial parameter of the stability 
of the A15 lattice and should not be too low nor too 
high. The A15 island in the phase diagram appears to 
be centered around k B T/e « 1.4 and its temperature 
range spans about 2k B T/e at most, the maximum being 
at S w 0.27. The A15 lattice is a delicate structure: for 
all 5 and T, the pure A15 structure occurs only within a 
rather narrow range of density centered at n rs 0.85 with 
width An 0.05. 

(ii) In the density-temperature plane, islands of sta- 
bility of the intermediate-density solid phases (i.e., ex- 
panded FCC, expanded BCC, A15, and condensed BCC 
lattice) are elongated along the temperature axis, imply- 
ing that over broad ranges of temperature, the phase se- 
quence does not change dramatically with density. This 
departs from the low-temperature extrapolations of the 
results of perturbative theoryO the stripe-like topology 
of the phase diagrams presented in Fig. [ll] and Fig. [l^ is 
free of critical and/or triple points as well as the corre- 
sponding isostructural transitions. This is most likely a 
consequence of the rather broad shoulder. 

(iii) For shoulders as broad as necessary to stabilize 
the A15 lattice, the phase diagram is characterized by 
relatively small islands of stability of the intermediate 
phases: two-phase coexistence dominates the phase dia- 
gram. This feature appears to be specific to convex inter- 
particle potentials, whereas in case of concave potentials 
the regions of coexistence are typically much narrower O 

While the stability of the A15 lattice is certainly not 
limited to HCSS potentials, our findings indicate that at 
least in some systems it could coexist with FCC, BCC, or 
perhaps another lattice over a rather wide density range. 
As far as the structural identification of the samples is 
concerned, this simple fact may have important conse- 
quences. In the case of phase coexistence, the interpre- 
tation of X-ray measurements is difficult. In practice, 
the space groups of the coexisting phases can only be 
determined unambiguously by varying the external pa- 
rameters and moving from one island of stability across 
the coexistence region to the other island of stability. 
Given the stripe-like topology of the phase diagram, the 
parameter to be varied should be the density and not 
temperature. Even with this proviso, the identification 
can be difficult because of a considerable overlap between 
the diffraction peaks of the structures involved. For ex- 
ample, the difference between the patterns of A15 and 
BCC lattice is not very striking - the reflections being 
at y/2, y/4, y/E, y/E, . . . and at yM, y/i, y/E, . . . for the A15 
and BCC lattices, respectively^ - and could be masked 
by the form factor of the particles,-which usually falls off 
quite rapidly with the wavevectorE3 These shortcomings 
could be overcome by complementing the X-ray studies 
by calorimetric measurements. 

Last but not least, we note that the relative rarity of 



the A15 lattice in real colloidal crystals may be caused by 
other, non-equilibrium mechanisms that can slow down 
the formation of a pure A15 colloidal crystal. For ex- 
ample, it is conceivable that the relatively simpler BCC 
lattice is kinetically favored over the A15 lattice in ex- 
periments which evaporate solvent to form crystallites. 
Given that the free energies of the different lattices are 
typically rather small, the equilibration can take very 
long and the actual ground state may be hard to-observe. 

Nevertheless, along-with recent experimental^ as well 
as theoretical studiesE2l our preliminary results point to 
the necessity of extending the phase diagram of soft 
spherical particles by the A15 lattice and possibly other 
non-close-packed lattices. At the same time, the more in- 
tricate phase diagram may be more difficult to determine 
and understand. 



IV. CONCLUSIONS 

In this study, we have extended the geometrical inter- 
pretation of the free energy of weakly interacting classi- 
cal particles, and we have complemented the well-known 
close-packing rule with a minimal-area rule. The incom- 
patibility of the two principles leads to frustration which 
gives rise to range of possible equilibrium structures, de- 
pending on the relative weight of the two terms. Our pro- 
posed scheme is a computationally simple, zeroth-order 
description of soft-sphere crystals which is more transpar- 
ent than detailed numerical models, most notably molec- 
ular modeling. 

As such, our theory provides a robust insight into the 
self-organization of such objects, which should be useful 
for the engineering of colloidal crystals. The relevance of 
these universal guidelines is as broad as the use of collejda 
themselves, ranging from photonic bandgap crystals£3'E£l 
to micro- or mesoporous materials-used for chemical mi- 
croreactors and molecular sieves. To meet the de- 
mands of a particular application, these designer mate- 
rials must be characterized by a given lattice constant, 
symmetry, and mechanical properties, and, in the case of 
porous structures, void size and connectivity. All these 
parameters can be controlled by tuning the structure and 
size of the (template) colloidal particles and the interac- 
tion between them, and our model establishes a semi- 
quantitative relationship between particle geometry and 
bulk material properties. 

We envision this work to be extended in several di- 
rections. One problem to be addressed is to locate the 
A15 lattice within the phase diagram using Monte Carlo 
analysis, starting with a square-shoulder potential but 
also employing less generic short-range potentials. In ad- 
dition, we will further explore the analogy between col- 
loidal crystals and soap froths in view of the geometrical 
approach that we have adopted here. Another interesting 
aspect of future work could be to use the model to de- 
rive some of the mechanical properties of colloidal crys- 
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tals, such as the shear and Young moduli. One could 
also study the stability of non-cubic lattices, which have 
been mostly neglected so far, in an attempt to understand 
non-spherical colloidal particles. Work along these lines 
should lead to easily verifiable predictions and a deeper 
insight into the physics of colloids. 
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